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Abstract. 

Both regular (the zero-momentum ghost dressing function not diverging), also named decoupling, and critical (diverging), 
also named scaling, Yang-Mills propagators solutions can be obtained by analyzing the low-momentum behaviour of the 
ghost propagator Dyson-Schwinger equation (DSE) in Landau gauge. The asymptotic expression obtained for the regular 
or decoupling ghost dressing function up to the order ff(q 2 ) fits pretty well the low-momentum ghost propagator obtained 
through the numerical integration of the coupled gluon and ghost DSE in the PT-BFM scheme. Furthermore, when the size of 
the coupling renormalized at some scale approaches some critical value, the PT-BFM results seems to tend to the the scaling 
solution as a limiting case. 
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INTRODUCTION 

The low-momentum behaviour of the Yang-Mills propa- 
gators derived either from the tower of Dyson-Schwinger 
equations (DSE) or from Lattice simulations in Landau 
gauge has been a very interesting and hot topic for the 
last few years. It seems by now well established that, 
if we assume in the vanishing momentum limit a ghost 
dressing function behaving as F(q 2 ) ~ (q 2 ) aF and a 
gluon propagator as A(q 2 ) ~ (q 2 ) aa ~ l (or, by following 
a notation commonly used, a gluon dressing function as 
G(q 2 ) = q 2 A(q 2 ) ~ {q 2 ) aa ), two classes of solutions may 
emerge (see, for instance, the discussion of refs. [1, 2]) 
from the DSE: (i) those, dubbed "decoupling", where 
a/7 = and the suppression of the ghost contribution 
to the gluon propagator DSE results in a massive gluon 
propagator (see [3, 4] and references therein); and (ii) 
those, dubbed "scaling", where ^ and the low- 
momentum behaviour of both gluon and ghost propaga- 
tors are related by the coupled system of DSE through 
the condition 2 Off + 0£<3 = implying that F 2 (q 2 )G(q 2 ) 
goes to a non-vanishing constant when q 2 — > (see [5, 6] 
and references therein). 

Lattice QCD results appear to support only the mas- 
sive gluon (ac = 1) or scaling solutions (see [7] and 
references therein), and also pinching technique results 
(see, for instance, [8, 14] and references therein), re- 
fined Gribov-Zwanziger formalism (see [9]) or other ap- 
proaches like the infrared mapping of A0 4 and Yang- 
Mills theories in ref. [10] or the massive extension of the 
Fadeev-Popov action in ref. [11] appear to point to. In 
the present note, we briefly review the work of ref. [13], 
which extended the previous studies of refs. [1,2, 12], by 



the analysis of the results obtained by solving the cou- 
pled system of Landau gauge ghost and gluon propaga- 
tors DSE within the framework of the pinching technique 
in the background field method [14] (PT-BFM) 



THE TWO KINDS OF SOLUTIONS OF 
THE GHOST PROPAGATOR 
DYSON-SCHWINGER EQUATION 

As was explained in detail in refs. [2, 12, 13], the 
low-momentum behavior for the Landau gauge ghost 
dressing function can be inferred from the analysis of 
the Dyson-Schwinger equation for the ghost propagator 
(GPDSE). That analysis is performed on a very general 
ground: one applies the MOM renormalization prescrip- 
tion, F R (p. 2 ) = ix 2 A fl (ix 2 ) = 1, where ji 2 is the subtrac- 
tion point, chooses for the ghost-gluon vertex, 

rf c (-q,k;q-k) = ig f bc {q v H x {q,k) 

+ {q-k) v H 2 (q 7 k)) (1) 

to apply this MOM prescription in Taylor kinemat- 
ics (i.e. with a vanishing incoming ghost momentum) 
and assumes the non-renormalizable bare ghost-gluon 
form factor, H\(q,k) = Hi, to be constant in the low- 
momentum regime for the incoming ghost. Then, the low 
momentum-behaviour of the ghost dressing function is 
supposed to be well described by 

F R {q 2 ) = A(n 2 )^ ' [ 1 + "-J . ( 2 ) 



and that of the gluon propagator by 
B(p 2 ) n( " 2 
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and this, after solving asymptotically the GPDSE, finally 
left us with: 
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if a/r = 0, where 



a T (0)=M : 
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-F 2 (0)A R (0) . 



(6) 



It should be understood that the subtraction momentum 
for all the renormalization quantities is ji 2 . The case 
(Xf 7^ leads to the so-called scaling solution, where the 
low-momentum behavior of the massive gluon propaga- 
tor forces the ghost dressing function to diverge at low- 
momentum through the scaling condition: 2a/r + = 
(a<j = 1 is the power exponent when dealing with a mas- 
sive gluon propagator). As this scaling condition is veri- 
fied, the perturbative strong coupling defined in this Tay- 
lor scheme [15], (Xt = gj/(4n), has to reach a constant 
at zero-momentum, 



Or(0) = 
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as can be obtained from Eqs.(3,5). The case a/r = cor- 
responds to the so-called decoupling solution, where the 
zero-momentum ghost dressing function reaches a non- 
zero finite value and eq. (5) provides us with the first 
asymptotic corrections to this leading constant. This sub- 
leading correction is controlled by the zero-momentum 
value of the coupling defined in eq. (6), which is an ex- 
tension of the non-perturbative effective charge defini- 
tion from the gluon propagator [16] to the Taylor ghost- 
gluon coupling [17]. 



The "critical" limit 

In ref. [13], the solutions of the coupled DSE system 
in the PT-BFM scheme (with H\ = 1 for the ghost-gluon 



vertex), numerically integrated for many values of the 
coupling at the renormalization point p 2 as a boundary 
condition, were studied. They were shown to behave 
asymptotically as eq. (4) predicts for the decoupling DSE 
solutions. There also appeared to be a critical value of 
the coupling, a^t = «(/i 2 ) ~ 0.182 with ji = 10 Gev, 
above which the coupled DSE system does not converge 
any longer to a solution. As a matter of fact, we know 
from refs.[2, 13] that the scaling solution implies for the 
coupling 



Ofcrit — 
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4k 2N c A 2 (ji 2 )B{jX 2 ) 
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where B{ji 2 ) and A(p 2 ) defined by Eqs. (3,4). This is 
also shown in ref. [1], where only the ghost propagator 
DSE with the kernel for the gluon loop integral is ob- 
tained from gluon propagator lattice data. In the analy- 
sis of ref. [1], a ghost dressing function solution diverg- 
ing at vanishing momentum appears to exist and verifies 
eqs. (4,8), while regular or decoupling solutions exist for 
any a < (Xcrit- In ref. [13], a more complete analysis is 
performed: first by studying the solutions for many dif- 
ferent values of the coupling, a = a(/i 2 ), of a coupled 
DSE system; and then by showing that the ghost dressing 
function at vanishing momentum, F(0,n 2 ), is described 
by the following power behaviour, 



F(0) ~ (a crit -a(M 2 )) 



(9) 



where k(/x 2 ) is a critical exponent (depending pre- 
summably on the renormalization point, fx 2 ), supposed 
to be positive and to govern the transition from decou- 
pling (a < Ofcrit) to the scaling (a = a cr i t ) solutions; 
and where we let Ofcrft be a free parameter to be fit- 
ted by requiring the best linear correlation for log[F(0)] 
in terms of log[a cr j t — a]. In doing so, the best corre- 
lation coefficient is 0.9997 for a cr ; t = 0.1822, which is 
pretty close to the critical value of the coupling above 
which the coupled DSE system does not converge any 
more, and k(h 2 ) = 0.0854(6). This can be seen in fig. 1, 
where the log-log plot of Fr(0) in terms of a cr ; t — a 
is shown and the linear behaviour with negative slope 
corresponding to the best correlation coefficient strik- 
ingly indicates a zero-momentum ghost propagator di- 
verging as a — > Ofcrit. Nevertheless, no critical or scaling 
solution appears for the coupled DSE system in the PT- 
BFM, although the decoupling solutions obtained for any 
a < a cr i t = 0.1822 seem to approach the behaviour of a 
scaling one when a — > a cr ; t . This is well understood in 
ref. [13], where the gluon propagators obtained from the 
coupled DSE system in PT-BFM were also found to obey 
the same critical behaviour pattern as the ghost propaga- 
tor, when approaching the critical value of the coupling. 
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FIGURE 1. Log-log plot of the zero-momentum values of 
the ghost dressing function, obtained by the numerical inte- 
gration of the coupled DSE system in the PT-BFM scheme, 
in terms of 0^ — a. a = = lOGeV), the value of the 
coupling at the renormalization momentum, is an initial con- 
dition for the integration; while a cr ; t is fixed to be 0. 1 822, as 
explained in the text, by requiring the best linear correlation. 

CONCLUSIONS 

The ghost propagator DSE, with the only assumption of 
taking H\ (q,k) from the ghost-gluon vertex in eq. (1) to 
be constant in the infrared domain of q, can be exploited 
to look into the low-momentum behaviour of the ghost 
propagator. The two classes of solutions named "decou- 
pling" and "scaling" can be indentified and shown to de- 
pend on whether the ghost dressing function achieves a 
finite non-zero constant ((fy? = 0) at vanishing momen- 
tum or not ((Xf ^ 0). The solutions appear to be dialed 
by the size of the coupling at the renormalization mo- 
mentum which plays the role of a boundary condition for 
the DSE integration. When applying a model with a mas- 
sive gluon propagator, the decoupling low-momentum 
behaviour of the ghost propagator results to be regu- 
lated by the gluon mass and the Taylor-scheme effec- 
tive charge at zero momentum and successfully describes 
the low-momentum ghost propagator computed trhough 
the numerical integration of the coupled gluon and ghost 
DSE in the PT-BFM scheme. The zero-momentum ghost 
dressing function is then shown to tends to diverge when 
the value of the coupling dialing the solutions approaches 
some critical value. Such a divergent behaviour at the 
critical coupling corresponds to a scaling solution where, 
if the gluon is massive, a/r = —1/2. 
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